Evidence for a low density Universe 
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X. 

^ , The motions of galaxies can be used to constrain the cosmological density param- 

eter, Q, and the clustering amplitude of matter on large scales. The mean relative 
velocity of galaxy pairs, estimated from the Mark III survey indicates Q = O.SSIq 25- 
If the clustering of galaxies is unbiased on large scales, = 0.35 it 0.15, so that an 
unbiased Einstein-de Sitter model (Q, = 1) is inconsistent with the data. 

The mean relative velocity for a pair of galaxies at positions fi and r2 is = Hf, where 
r = f*! — and the constant of proportionality H = 100 h km s^^Mpc^^ is the Hubble parameter 
(1,2). The quantity 0.6 < h < 1 parametrizes uncertainties in H measurements. This law is an 



idealization, followed by real galaxies only on sufficiently large scales, corresponding to a uniform 
mass distribution. On smaller scales, the gravitational field induced by galaxy clusters and voids 
generates local deviations from the Hubble flow, called peculiar velocities. Correcting for this effect 
gives = Hr + Vi2 r/r. The quantity Ui2(r) is called the mean pairwise streaming velocity. In the 
limit of large r, = 0. In the opposite limit of small separations, 1X12 (r) = (virial equilibrium). 
Hence, at intermediate separations fi2 < and we can expect to observe gravitational infall, or 
the "mean tendency of well-separated galaxies to approach each other" (3). In a recent paper we 
derived an expression, relating V12 to cosmological parameters in another, using Monte Carlo 
simulations we showed how can be measured from velocity-distance surveys of galaxies (5). 
Our purpose here is to estimate t'i2(r') from observations and constrain the cosmological density 
parameter, 

The statistic we consider was first introduced in the context of the Bogolyubov-Born-Green- 
Kirkwood (BBGKY) kinetic theory describing the dynamical evolution of a self-gravitating collection 
of particles ( 3, 6). One of the BBGKY equations is the so called pair conservation equation, relating 
the time evolution of V12 to ^(r) — the two-point correlation function of spatial fluctuations in the 
fractional matter density contrast (3). Its solution is well approximated by (4) 

v,2{r) = -ii7rJ^0-6f(r)[l + a|(r)], (1) 

= _ 3 j;; ^jx) dx 
^^"^ - r3[l + e(r)] ' 

where a = 1.2 — 0.667, 7 = —{dhi^/dhir)^=i and f2 is the present density of nonrelativistic 
particles. The equations above have been obtained by interpolating between a second-order per- 
turbative solution for 1^12 (r) and the nonlinear stable clustering solution. For a particle pair at 
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separation r, the streaming velocity is given by 

Vi2{r) = {{V1-V2) ■f)p = {{vi - V2) ■ fwi2) , (3) 

where w;i2 = {1 + 6i){l + 62) [1 + ^(r)]~^ is the pair-density weighting, Va and 5a are the peculiar 
velocity and fractional density contrast of matter at position Va, A = 1,2..., the separation 
r = |ri — is fixed for all pairs, the hats denote unit vectors, and ^(r) = {6162) ■ The expression in 
square brackets in the definition of W12 ensures that {W12) = 1 and the pairwise velocity probability 
density integrates to unity. Note that the pair-weighted average, (• • differs from simple spatial 
averaging, (■■■), by the weighting factor Wi2- The pair- weighting makes the average different from 
zero, unlike the volume average {vi — V2) = 0, which vanishes because of isotropy. 

Our approximate solution of the pair conservation equation was successfully tested against 
N-body simulations in the dynamical range ^ < 10^ (4>V- I* valid for universes filled with non- 
relativistic particles and it is insensitive to the value of the cosmological constant (2,4)- Eqn.(l) 
was derived under the additional assumption that the probability distribution of the initial, small- 
amplitude density fluctuations was Gaussian. 

Until now we have also implicitly assumed that (i) the spatial distribution of galaxies traces 
the mass distribution and that (ii) Vi2{r) for the galaxies is the same as for the matter. If the 
galaxies are more clustered than mass, condition (i) is broken and we have "clustering bias" . The 
galaxy two-point correlation function is close to a power law, ^^^'(r) oc r~'^, over three orders of 
magnitude in separation r (8). This is not true for the mass correlation function ^(r) in structure 
formation models of the cold dark matter (CDM) family (7). To reconcile theory with observation, 
one has to introduce a measure of bias that depends on separation and cosmological time, t: 
b'^{r,t) = ^^^\r,t)/^{r,t). Because of the pair-density weighting, clustering bias can in principle 
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induce "velocity bias" in a way similar to systematic error propagation. This is certainly true in 
the most simplistic of all biasing prescriptions - the "linear biasing" , under which 6 is a constant 
and, moreover, 5^^^ = b6. The expression for can be obtained from Eqn.(^) by formally 
replacing the weighting function 1012(61,62) with ti'i2(5P\ In the linear limit, ^ ^ 1, we get 

Vi2^ = bvi2 (9), in qualitative agreement with recent N-body simulations, which considered a whole 
range biasing prescriptions, allowing nonlinear and/or non-local mapping of the mass density field 
onto 6^^^^ (10). However, there are also simulations which show exactly the opposite: although the 
galaxies do not trace the spatial distribution of mass, pairs of galaxies behave like pairs of test 
particles moving in the gravitational field of the true mass distribution, and vf2'^{r) = fi2(r) (11)- 
Direct measurements of vff'{r) can help us decide which simulations and biasing schemes are more 
believable than others. Indeed, one can measure vff' for different morphological classes of galaxies. 
The linear bias model predicts Vi2^^V'^i2^^^ = h^^^ /b^^\ where the superscripts refer to elliptical (E) 
and spiral (S) galaxies. Observations suggest ^^^^/ft*-^-* ~ 2 and h^^^ ~ 1 (13). Hence, one expects 
Ui2^^V^i2^^'' ~ 2 if the linear bias model is correct and 1^12 '•^V^ia^^'* = 1 in the absence of velocity 
bias. 

Measurements of t'i2(?') can be also used to determine Vl. Indeed, if the mass correlation 
function is well approximated by a power law, ^(r) oc r"'*', V12 at a fixed separation can be 
expressed in terms of and the standard normalization parameter - Ug. The latter quantity is the 
root-mean-square contrast in the mass found within a randomly placed sphere of radius 8 h~^Mpc. 
Unlike the conventional linear perturbative expression for 1^12 (?") oc ^^'^ r'^~^ , our nonlinear 
Ansatz provides the possibility of separating from by measuring V12 at different values of r 
[see the lowermost panel in Fig.l below; see also ref. (12) ]. 
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We will now describe our measurements. The mean difference between radial velocities of a 
pair of galaxies is {sa — Sb )p = ^12 ^ • (^a + ^b)/2, where Sa = Ta • Va and r = Va — fs- Here as 
before, the latin subscripts number the galaxies in the survey, A,B = 1,2 .. . To estimate V12, we 
minimize the quantity X^(fi2) = E^,s [{sa - Sg) - Pab ^^12/2]^ , where Pab = f • (f^ + rs) and 
the sum is over all pairs at fixed separation r = |f^4 — rgl- The resulting statistic is (5) 

2J2{sa- Sb)pab 

= — — ■ 

Monte-Carlo simulations show that this estimator is insensitive to biases in the way galaxies are 
selected from the sky and can be corrected for biases due to errors in the estimates of the radial 
distances to the galaxies (5). The survey used here is the Mark III standardized catalogue of galaxy 
peculiar velocities (14,15,16). It contains 2437 spiral galaxies with Tully-Fisher (TF) distance 
estimates and 544 ellipticals with D„ — a distances. The total survey depth is over 120 h~^Mpc, 
with homogenous sky coverage up to 30 h-^Mpc. The inverse TF and IRAS density field corrections 
for inhomogeneous Malmquist bias in the spiral sample agree with each other and give similar 
streaming velocities, with lognormal distance errors of order aind ~ 23%. For the elliptical sample, 
cind ~ 21% and the distances assume a smooth Malmquist bias correction (17). 

The estimates from the spiral and elliptical are remarkably consistent with each other (Fig.l), 
unlike previous comparisons using the velocity correlation tensor (18,19). For a velocity ratio 
R = Vi2^^yvi2^^^ = 1, we obtain ~ 1, while for i? = 2 the = 2.1. The most straightforward 
interpretation of this result is that there is no velocity bias and the linear clustering bias model 
should be rejected. Its static character and the resulting failure to describe particle motion, induced 
by gravitational instability was pointed out earlier on theoretical grounds (20). Our results can, 
however, be reconciled with linear bias model if it is generalized to allow scale-dependence, b = b{r). 



Biasing factors for both galaxy types can be arbitrarily large at small separations, where ^(r) ^ 1, 
if biasing is suppressed at large separations, where |^(r)| < 1. Indeed, in the nonlinear limit 

■Wi2(6(5i, Ma) h'^5i52/h'^i = 5i52/i, and hence vff{r) Ui2(r). 

We obtained an estimate of cjg and $7 from the shape of the v-i2{t) profile as follows. We assumed 
that the shape of the mass correlation function ^(r) (but not necessarily the amplitude) is similar 
the shape of the galaxy correlation function estimated from the APM catalogue (8), consistent 
with a power law index 7 = 1.75 it 0.1 (the errors we quote are conservative) for separations 
r < 10 h-^Mpc. Given the depth of the Mark III catalo gue we expect the covariance between 
estimates of v-i2{r) to be only weakly correlated at r < 10 h^^Mpc; we use N-body simulations to 
determine the covariance of the estimates over this range of scales and use a minimization to 
obtain the 1-a constraints: Cg > 0.7 and = 0.35^q'25. Fixing = 1 we obtain VL = 0.35 it 0.15 
(Fig. 2). 

We can obtain a more conservative constraint on cjg and $7 by examining a v-12 at a single 
separation, r = r* = 10h~^Mpc. Substituting r = and ^(r) oc r"^-^^ into Eqns.(|l|)-(P), we get 

v^2{r*) = - 605 (jg^OO-S (1 + 0.43CJ82) /(I + 0.380-82)2 km/s . (5) 

The above relation shows that at r = r*, v-12 is almost entirely determined by the values of two 
parameters: Ug and Vi. The uncertainties in the observed 7 lead to an error in eq. (^) of less than 
10% for (jg < 1. In fact, at this level of accuracy and at this particular scale, our constraints 
depend only on the value of and the overall normalization cjg but do not depend on other model 
parameters, such as the shape of ^(r). The streaming velocity, v-^2{f*) depends on ^(r) only at 
r < r=K, so unlike bulk flows, it is unaffected by the behavior of ^(r) at r > [compare our Eqn. (1) 
with Eqn. (21.76) in ref. (2)\. Moreover, the dominant contribution to Wi2(?'*) comes from ^{r^) - 
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an average of ^(r) over a ball of radius r*, so the details of the true shape of ^(r) at r < r^, have 
little effect on 1712 (r*) as long as cjg (and hence, the volume-averaged ^) is fixed. Hence, Eqn.(|5|) 
can provide robust limits on cjg and even if the assumption about the proportionality of ^(r) 
to the APM correlation function is dropped. This statement can be directly tested by comparing 
predictions of Eqn.(^) with predictions od CDM-like models, all of which fail to reproduce the pure 
power-law behavior of the observed galaxy correlation function. When this test was applied to four 
models, recently simulated by the Virgo Consortium (1)^ we found that for fixed values of cjg and 
the predictions based on Eqn.(P) were within < 6% of the Vi2(?"*), obtained from the simulations 

m 

The measured value, — Ui2(?'*) = 28OI53 km/s (Fig.l), is inconsistent with = 1 and = 1 
at the 99% confidence level. 

Our results are compatible with a number of earlier dynamical estimates of the parameter 
/? = rj^'^cTg [/3 is sometimes defined as fh^ but Cg ~ 1/6 and the two definitions differ at most 
at the 10% level, see, for example ref. ^.^j ]. A technique, based on the action principle (22) gives 
/? = 0.34 lb 0.13; comparisons of peculiar velocity fields with redshift surveys based on the integral 
form of the continuity equation (called velocity-velocity comparisons) typically give /? = 0.5 — 0.6 
(23-26). Of the velocity-velocity comparisons, the one with the smallest error bars is the VELMOD 
estimate: /3 = 0.5 it 0.05 (24)- This constraint has several advantages over others; in particular, 
it correctly takes into account cross calibration errors between different Mark III subcatalogues. 
To illustrate the consistency of our results with velocity- velocity studies we will now compare our 
limits on c7g and fl, derived from the shape of Vi2{r) for a range of separations with constraints from 
our measurement of i'i2(10h^^ Mpc) alone, combined with limits from VELMOD (Fig. 2). Again we 
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find that a low-O universe is favoured: U < 0.65 and erg > 0.7. The concordance region overlaps 
with the constraint derived from our measurements of fi2(r). 

Our results disagree with the IRAS-POTENT estimate, /? = 0.89 ± 0.12 (27). The IRAS- 
POTENT analysis is based on the continuity equation in its differential form; it uses a rather 
complicated reconstruction technique to recover the full velocity field from its radial component. 
The reason for the disagreement is not clear at present. One can think of at least two possible 
sources of systematic errors in the IRAS-POTENT analysis: (i) the reconstruction scheme itself (e.g., 
taking spatial derivatives of noisy data), and (ii) the nonlinear corrections adopted. The nonlinear 
corrections diverge like in the limit $7^0 (27). By contrast, the accuracy of the nonlinear 

corrections for the velocity- velocity is insensitive to CI. The velocity- velocity approach is also simpler 
than IRAS-POTENT because it does not involve the reconstruction of the full velocity vector from its 
radial component measurements (although both approaches do require a reconstruction of galaxy 
positions from their redshifts). Note that our method is direct, not inverse: it does not require any 
reconstruction at all. 

Finally, there is a potential caveat in the "no velocity bias" assumption in our own analysis. 
Although this assumption is based on empirical evidence from the two sets of galaxy types, the 
observational data is noisy and involves non-trivial corrections for Malmquist bias which could 
affect the two samples differently. Application of our approach to different data sets will clarify 
these issues. If contrary to our preliminary results, the streaming velocity turns out to be subjected 
to bias after all, such a finding may affect our estimates of erg and fl, based on the shape of the Vi2{r) 
profile but not our rejection of the unbiased Einstein-de Sitter model. In this sense our differences 
with the IRAS-POTENT analysis do not depend on the presence or absence of velocity bias. 
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The advantages of the new statistic we have used here can be summarized as follows. First, 
Vi2 can be estimated directly from velocity-distance surveys, without subjecting the observational 
data to multiple operations of spatial smoothing, integration and differentiation, used in various 
reconstruction schemes. Second, unlike cosmological parameter estimators based on the acoustic 
peaks, expected to appear in the cosmic microwave background power spectrum (28), the O esti- 
mate based on Vi2 is model- independent. Finally, our approach offers the possibility to break the 
degeneracy between Q, and cjg by measuring the Vi2{r) at different separations. 
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Figure 1: The streaming velocities of 2437 spiral galaxies (top panel) and 544 elliptical galaxies (center 
panel) estimated from the Mark III catalogue. The error bars are the estimated la uncertainties in the 
measurement due to lognormal distance errors, sparse sampling (shot noise) and finite volume of the sample 
(sample variance). The error bars were estimated from mock catalogues described in ref. (5). The small 
sample volume also introduces correlations between measurements of Vi2{r) at different values of r. To guide 
the eye, and show that although the two samples have different noise levels (because of much smaller number 
of galaxies in the elliptical sample), the Vi2{r) signal in both cases is similar, we also plot Vi2{r) calculated 
from equation (1) for a ^ oc r"^ ''^ power-law model with cTg = 1.25 and fl = 0.3. Three theoretical Vi2{r) 
curves are plotted (bottom panel) with ^ oc r~^'^^, CTgfi^-^ = 0.7 and cTg = 0.5 (solid line), 1 (dotted line) 
and 1.5 (dashed line). These curves show how measurements of Wi2(r) can break the degeneracy between O 
and CTg. 

Figure 2: The blue region constrains the viable values of O, the fractional mass density of the universe, 
and cTg, the variance of mass fluctuations at r = 8 h~^Mpc, from the combination of the constraints on the 
streaming velocities (red region) and P = (green region). The streaming velocities are constrained at 

r = 10/i~^Mpc from the Mark III catalogue of peculiar velocities and /? is measured using the VELMOD 
comparison between the Mark III catalogue and the velocity field inferred from the IRAS redshift survey. 
The dashed line defines the 1-a obtained from comparing expressions 1 and 2 with the Mark III catalogue 
from 2h-iMpc to lOh'^ Mpc. 
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